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:/ffiSTIlAGT 


AlgoiithEis for cioss hatched aa::d half toned representation 
of surfaces has heen developed i Cross hatched, representation 
algorithn is applicable to surfaces obtained by continuous linear 
transfornation of, function of one variable. For this family of 
surfaces our algorithm is much faster than the already existing 
algorithms. For eliminating the hidden part of the surface vra 
have modified the normally used ¥atkin*s masking algorithm* Half- 
toned representation algorithm is applicable to surfaces of the form 
Z = f(x,y). Here the rectangle in which the surface is defined is 
divided into subreotangles depending on the partial derivatives 
of z at left base corner of the rectangle then the points on the 
centra of the rectangles are projected on the view plane, to get the 
shadded picture of the surface. 
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1 . IITTROBUCTIOSr 


This is a well Icnovm sa 3 /ing that "A picture is vrorth a thousand 
words" 5 that this proverb holds today in the fields of physical sciences 
and engineering is demonstrated in the abundance of illustrations and 
figures contained in technical reports. 

The two common methods of representing a surface is contour 

plots and planar projections. Contour plots depicts the qualitative 
properties of the surfaces and are useful) if' numerical information 
is to 'be extracted from the figure. A planar projection depicts the 
surface in such a manner that its properties are most easily visualized. 
Further, planar projection has two common categories, cross hatched 
representation and continuous shedding or what is called a half tone 
reprosenAtion. In cross hatched representation mutua?.ly intersaoting 
family of contours (normally one with fixed X-values, and second with 
fixed 1— values) are drawn to give the illusion of three dimensionality, 
whereas in half-tone representation points witn varying density sire 
plotted to give the 5-1 view of the surface. Both the categories of 
projections have two kinds, orthographic projection where lines 
perpend iqular to the viev/ plane from the points on the surface are 
drawn which corresponds to the observer at infinite and perspective 
projection where line from the points on surface is joined to the point 

of observation, and in both the oases points or lines are drawn on the 
view plane where these set of lines intersects the plane. To give 
a proper 5-1 view of the surface, those parts of it ?;hich are hidden 
from the observer are removed while drawing 'the surface. 
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We have devQlopad algoiiithniE for cross hetched ai'd half tong repr©- 
sjaataiion- ' surfaces with orthographic projection. As there are 
already many algorithms existings there oan not be any justification 
of giving a new algoritlrm v/ithout any si?:nificant improvement in the 
time or spans oompoxity of them. The time complexity of the existing 

algorithms for cross hatched representation is dependent on both the 
complexity of the surfaces as well as the number of points considered 
on them, because if the surface is 2 = f(x,y) the function f is 
evaluated for different values of x and y. Our algorithm is applicable 
to only a particular family of surfaces which are obtained using some 
continuous transformations on functions of one variable, for which- we 
do away with the need of calculating (x,y,z) for different points chosen 
on the surface. We calculate (x,y,z) for points on only a particular 
ourwe lying on tho surface and remaining points on the surface are 
obtained from these points using the transformation. This makes the 
time complexity of our algorithm independent of the complexity of the 
surface and depends only on the number of points taken on it. More- 
over, if the equation of the surface is f(x,y,z) = 0, where z is not 
easily separable, then the existing algorithms require to solve the 
equation f(x,y,z) = 0 for fixed va.lues of X and Y, making it even more 
complicated to determine (x,y,z), whereas our algorithm even though 

not applicable to any general f(x,y,z) = 0, generates many surfaces 
using transformations which can not be explicitly solved for z. (For 
example, see the surfaces in Figures 5 and 6.) The time complexity 
of our algorithm is dealt in Chapter 5* For the surfaces where our 
algorithm is applicable it is much faster than the existing algorithms. 



3 


I.' or Gliminating hidden lines v/e have modified the normally used 
Yfatkin's algorithr to make it faster and applicable to discontinuous 
curves . 



2. Et/iSKIl'IG ALGCRITBI'IS (Hidden-line Elimiiiation.) 

In genexal for a given surface and angle of vie-.; not all parts 
of the surface are visible, to gpt a prox)er vie-? it is required to 
eliminate that part -v/hich is hidden from the observer. One of the 
commonly used algorithms for eliminating the hidden part is by Yfetkins, 
vrhich takes a curve on the surface as input, eliminates the hidden 
parts of it and draws the rest. Our algorithm is a modification of 
it. Watkin's algorithm is briefly given bolov;, for detail see 
Referenoo 5 • 

2-1 . S ALGORITHM 

'Iho surface is scanned from front to rear* After proper rota- 
tions to ^t the required vieiv it is projected on the view plane 
(place perpendicular to the line of viev; vfhich in the case of ortho- 
graphic view is an unique plane). A proper coordinate system is 
chosen on the viav; plane. Two arrays ma-sk loi.- and mask high a-ro 
maintained , 5 which respectively contains for each possible pen position 
across the page the minimum and maxim-mn values attained so far, 
for the points projected on the vie-;; plane. Hask lev/ and mask high 
define two curves. That part of the segment which lies botvreen these 
two curves would bo hidden. This algorithm to sec which part of the 
segment lies between the two masking curves, determines for each 
possible pen position across the 'page on the segment, the Y-values 
and is oompared with that of corresponding mask high and mask low, 
if it lies between these values then this point ’.Tould be hidden. 

And the pen is moved to this now point without drav;ing any line. 
Otherwise if it is not hidden then a straight line is dravm from the 


4 
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pi’QVi^ Jis pen position to this point and tho masking ari^ays suce updated, 
iis the check is m.'^de at each and eveiy pen position an i acooidingly 
lino drawing loutine is called to dnaw oithex a olank line or solid 
line, time required for masking a.s well as dra.vir.g the lime is consi- 
derably large, consequently making the surface drawing rather slow. 

2-2. MOlIFICATIQg OF WJKUVS ALGORITHM 

We can do aiTay with tho need of checking for hidden point for 
faoh pen position with only a slight effect on the accuracy of hidden 
line elimination- 1 ■' '■ - 

Let the line segment- be given in an interval Tajh’!, and the 
coordinates of the end points of tho line segment be and (bjYg) 

respectively. The values of mask lev; and mask high at a and h are 
klLa, ICLb, MEa and IvIIIb respectively. Approximating that mask low and 
mask high vary linearly in the interval ia,H , the intersection of 
line segment [(a. (h^yg)] with i(a,l'Ila), (b,irLb)j' and 


|(a,I'iEa), (b,MHb)j is given holovfs 

(1 ) Intersection with mask lov7’, 

I Z1 = (itta * ) / (l-ILa - J£Lh + Y, - Y^ ) (1 ) 

I Y1 = a + (a-b) * I X- 1 / (y^ - Yg) (2) 

where (l X1 , IY1 ) is the intersection point. 

( 2 ) Intersection with ma,sk high - 

I X2 = (Ma Yg “ Jmb Y^ ) / (Mia - I.5I-Ib + Y^ - Y^ ) ( 5 ) 

I Y2 = a + (a-b) * I Xg / (Y^ - Y^) (4) 

where (lX2, IY2) is the intersection point. 


The above equations holds when the following conditions ans 


true 



Ca) l,!La - BILb + ^ 0 

(b) MHa - Ifflb + Y2 ~ / 0 

(c) Y^ - Y^ 0 

wTaen the above conditions hold then depending on the follouing 
co,s6s paxt or whole of the line segment would be visible 

(1 ) If a < I Z1 /..b and a <. I X2 Y.. b then the line 

segment intersects v/ith both the masking arrays and its 
portion lying betvraen I Z1 to I X2 would be hidden. 

(2. a) If a < I XI < b end (l X2 or I X2 > b) then the 

line segment intersects with only mask low. In that 
case if MLa ^ Y1 then the portion betmen I Z1 to b vrould 
be hidden else the portion between a to I XI would bo 
hidden. 

(2.b) If a<.I X2 < b and (l X1 <. a or I X1 >b) then the line 

segment intersects with only mask higb and if Ma Y1 then 
the portion hstv/eon 1X2 to b would bo 'hidden else the 
portion between a to I X2 would bo hidden. 

(5) If (iCLa <Y1 < ma and Iffib < 12 <, ICEIb) or 

(Y1 > lIHa and 12 > ffilb) or (yI <, 14La aaid Y2 < lILb) 

then the whole line segment would be visible. 

Cases \7hen one or more of conditions (a), (b) and (c) are false; 

(1 ) \’3hen only (a) is false - 

The line segment |(a,y^), (bjyg)! would be parallel 
to the mask lo’w segment. If y^ <C^ IvILa then the v/hole segment 

TiTould be visible, else if a ^I X2 <C, b then the analysis runs 

exactly as in the case 2.b above otherv/ise if Tl IIHa then 
the whole segment would be hidden. 

(2) l/djen only (b) is false, the analysis runs similar to the 


above case 
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(3) Wlnen (c) is false - 

line segment !^a,y^ ) , (b^y^)! would be parallel to 

the X-axis, I X1 sm& I X2 vrould Ise dstorinined by the 
Equation (1 ) and (5) respectively. I Y1 and I 12 vrould be 
equa.! to Y1 . And the hidden line analysis vrould be same 
as the case when all the three conditions hold. 

(4) Wisxi (a) and (b) are false - 

If then the whole line segment vrould. 

be hidden otherwise it vrould bo visible. 

These are the only four possible independent casos. Eamaining 
combinations are covered in the cases discussed above. After drav/ing 
each segment of line the masking arrays are appropriately updated. 

The above scheme for hidden line elimination works quite well 
when the masking curves are not very fast varying functions, otherwise 
also by choosing the interval |^,bl smaller it can be made applicable. 

The other problem v;ith Catkin’s algorithm is that it works only 
for continuous curves. Suppose that the partial function f is defined 
in the interval [X1 , X2j , end is discontinuous or undefined in 
certain subintervals. If it is passed to ¥atkin*s algorithm to draw, 
after eliminations the hidden part it vrould draw straight lines in the 
intervals whais it is undefined. Eor example in drawing a hyperboloid 
it would connect its vertices by a staight line, 'fhe proolem is 
solved very easily by first compbfcing the partial function f, by 
assigning a largo value ^4. all the points whore itwas undefined. 

And the masking algorithm for each given x in the interval 
checks whether f (x) is equal to b; , if it is equal to "’X. ■*'hsn 
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no line is f'xoiii ihG pisvious poiiit to this point. If f(x) 

is not equal to but at the previous point the value of the function 
f was then a blank line is dravm fron the last point v;hcre f was 
not ^'] to the present point. See figure 4 (a hyperboloid intersobting 
with a cylinder). 

2-3* BBLill lG fJ G CaiPLjXI ir Of 2J0 HiSKIHG ilGOllIfEIS 

Suppose the function f is continuous on the closed interval 
and is s^pproxiciated by M-1 line segaents and the nuiiiber of possible 
pen positions in the interval j is Im (N is noriaailsf much 

larger than II) . 

(1 ) V/atkin’s algorithn calculates the function f at all the 
possible pen positions II by linear interpola.tion and maices a chock for 
hidden point. Let the tine required to cheok ’.rhethcr a point is hidden 
or not be > and to calculate the value of the function f a-t a, point 
by linear interpolation bgr jj^ then the total tine required for masking 
by Watkin’s algorithm would be - 

I * + (it - II) * p" (5 ) 

normally M is of the order of 100 and IT afoout 1000. So the time 
complexity of ifetkin's algorithm is 0(ll) which is independent of 
the number of line segment the curve is made of, Y/hich means that even 
a constant function (f (x) = c) v^ill take as much tiiaa as any coriiplica4ed 
function would take. 

Moreover if making a call to draw a straight line takes TT tiiaa 
(which for small line segment is independent of its length), then tho 
the tine required to draw the curve would be TT * H which is aga.in 
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(2) In oui al^jorithn (nodification of Watkin’s algorithm) inter- 
section of each line segment is dotorminad with hoth the masking amays. 
If finding the intersection point and checking which part of the segment 
is above ' the masking function end v/hioh he low, takes on average &■ tine 
then finding which part ox the curve is hidden will re (pairs 2(il-1 ) * 9 
time, (ihe fa.ctor 2 i.s hocause -jq find the intersection point with 
hoth tho masking arrays). Moreover in updating the maskiig arrays we 
have to interpolate the function atmost at (N-M) pen positions. As 
each interpolation time is so tho total tine required hy our algorithm 
is - 

2(M-1 ) * 9 -1- (lI-M) (6) 

is usually much smaller than 9 and V . 

Comparing Equations (5) and (6) we see that second torn is common 
in hoth of than hut the first tarn in Vatkin’s algorithm contains H 
¥:hcre as our*s M. As said before M is much less than IT (iT is at least 
10 times of M), also 9 and Y aro of the same order, so the second 
algorithm is several tines faster than the first. Finally the time 
required to draw the curve in our algorithm is * Of“l)> which 
again is 0 (m) and much less than '.fetkin’s algorithm. 



3. SURPilCE DISPLAY iLGOfiimiS 


In the folloTfing sections tv/o different surface display algorithms 

are doscrihed. The first one diavz-s cross hatched surfaces, i7hich are 

obtained by continuous transformation of a function of single variabdn • 

The second one dravra half tone picture of the surface of form z = f(x,y). 

5-1 • DISPIAY OP A SUBPAGE OBTAIIED BY GSRBllALIZED OP A ODRYS 

ABOUT A LPE 

Generalised rotation as the name suggests is ge-neralisation of the 
ordinary rotation* In rotation the distance of a point from the line 
of rotation remains constant, nheroas in generalised rotation it varies 
T/ith the angle of rotation* To be precise the generalised rotation 

about a line and by an angle A is defined as follows s 

A 

let IB be the set of real numbers, ROT, - the ordinary rotation 

transformations about T and by an angle A* Then the generalised 

5 

rotation G. is a linear transformation on IB - 
A, < 

• 5 3 

such that 

W (X) 

where X €. IR^ and f is a continuous function from D, a compaot 
connected subset of R to IR such that f (o) = 1 • 
iis can be easily seen ordinary rotation is a particular cetse of 
the generalised rotation when f(A) = 1 ¥ afcIR. If a function y » H(x) 
is rotated about the X-axis using the generalised rotation function f 
then it can be easily -verified that the equation of the surface 
generated would be - 
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= (h(x)^ * (fOb:otan(z/y)))^ (1 ) 

It is intGiestinp io see that a plane suiiaoe can be obtained by the 
notation of a line using the gene2:0.1isGd notsiion function - 

f(i) = Cosec (A+P) c 
Assuming that the oquation of the line is 

Y = aX + b ox H(x) = ax + b 

and using Equation (l ) to denive the equation of the surf ace generalised s 

■f + Z = (a. X + b) * c Cosec (/occtan z/y +P) 
or 

(y^+Z^) Siji^(Acctan(z/y)+P) = (a X+b)^ * 
oi: 

(y^+Z^) * (Sin Ar:otan(z/y) * Cos P + Cos iicetan (z/y) * 

Sin P)^ = (a X + b)^ * 
or 

(y^ + Z^) * (Cos P * Z /-/yZ+^Z + Son P * y / 

= (ax + b)^ * 0 ^ 


or 

z Cos P + Y Sin P = (ax + b) * c 
vbich is G quoit ion of a plane '• 

The acLvantage of surfaces obtained by such generalised rotation is 
that there is an elegant iray of cross hatched representation using two 
families of functions, where the calculation of (x,y,z) for the points 


on then is very easy and. independent of the complexity of the surface' 
One of the family of functions is H(x) successively rotated (in thc- 


generalised sense) by an angle SA and peonseted on the’XfY plane ' ' 


63019i 


*». N». rxrs 






FjeUFE 1 : FC 0) *=1 , NO. OF POINTS-43Z0 
CPU TIME-12.69 SEC. 


FIGURE 2: FC 6) =1 . 0 

SURFACE OBTAINED BY THE ROTATION OF A 
SEMI CIRCLE. (SIMILER TO CYCLE TYRE) 
NK). OF POINTS-1216, CPU TiriE-»3.91 SEC. 


FIGURE 3; Ft 01 -EXF< -0^6025 _ ^ ^ . 

NO. OF F^INTS-2300. CPU TlhE>=6.B7 SEC 




NO?^ POINtI^KPi^* CYLINDER 

tw. Ut- PUINTS«B200. CPU TIME«18,26 SEC. 



FIGURE 6: FC©i«»1.0. H'YPEF?BOLjOID AND CYLirCER 
NO. OF POINTS-9427. CPU TiriE»21 ,09 SEC. 
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(assiMing that the ohseiiver is in the direction of Z~aocis at infinity 
in case the d inaction of obsenvation is diffenent from z then the sur- 


facQ is rotated to bring the Z-axis in the direction of obsortation), 


tmtil it has bc-on rotated by an angle of 560* The points on tho un- 
rotatod function is calculatad using Y = H(x) (z is zero for than), 
and tho generalised rotation matrix is determinGd, v/hich is 

successively applied to the function previously calculated, to give 
other members of the family. The second family consists of functions 
in planes parallel to the Y-Z plane lying at equal intervals in the 
domain of the function H(x), obtained by multiplying the- generalised 
rotation function f(A) by H(x). Hero again f(A) is determined, and 
all the functions of this family aro obtained by multiplying H(x) 

(x is difforont for difforont membors of tho family) to tho alroady 
calculated f(A). (in this caso the transformation is one dimonsional) • 
Tho gonoral mombor of tho first family is givon by 

ii 0 0 1 r X. I 

f I 1 

I 111 

P 1 ^(i) - f(A) * jo CosA SinA j )( j H(x^) j (2) 

I 0 -SinA Cosiy f 0 j 

■fch 

whoro F1j^(i) is tho coordinate of tho i point of that raonbor of tho 
family obtained by rotating H by angle A« x^*s and H(_x£)*s arc sane 
for all PI 

ii 

Tho -mombor 3 of tho second family arc obto,inGd by 


r 

J 


X 


P2^(9) = 


H(x)*f( 0 ) Gos e 
/ / 

[_H(x)*f(e) Sin ©J 


( 3 ) 



xnGSG itinctions arc passed to the navsking routino whore after 
proper rotation and oliiaination of tho hidden part the functions oro 
drawn. 

The flowchart of the a]gorithn is given 'bGlo^v’s 


Stop 

( 1 ) 

( 2 ) 


( 5 ) 


(4) 


STOP ).- 
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'Jhc coraplatc piograE is given in Appond ix-*1 . hrvG plotted 6 

suxfacQs using this algoiithra. The total nuuber of poiiri:s takon on the 

suifacG and the CPU tiiao is given v/ith 'inch figuno* Piguzc 1 oind Pigupe 

oix'clc 

2 aiG obtained by rotating a seni / about Z-cccia with f (xi) = 1 . In 
the first oaso the centre of tho cirolo lies on the Z-ojiis vrhereas in 
the second it is at (O.l , 2.0)^ tm.3.3 gEinerpiing a sphere .and a. .surfo-co 

similcu: to that of a cycle tyro respectively. Surfo.ee 5 ai:id 4 
obtained by rotating a straight line usiiig the gencralisod rotation 
function ]ilxp(-A/600) and (l+Sin 2A/2) respectively and the equation 
of surfaces 5 ^nd 4 s-rs given by:; 

Z = boolean (-2 Ln((y^ + Z^)/2.5)) and 

(Y^ + Z^)^ = 2.3(Y^+Z^ - IBS(YZ))^ respictively. 

'These are derived using Equation 1 • As we see the above functions 
cannot be explicitely solved for Z, so if we try to plot them using 
conventional method the above two equations will have io bo solved 
for each value of Y (fortunately they are indopendant of X). Biit our 
algoritto hojidly requires any computation, and draws then in about 7 
seconds. 

Analysis of the tiae complexity of tho algoritlra 

Consult the flowchart on Page 15 and the program in Append.ix 1. 
Suppose wa consider IT points on the function H(x) and the time, roquirad, 
to compute tl(x) is t1 . Then step one to.kcs IT*t1 time. In step tvro 
the U points calculated above arc rotated successively by an angle of 
''i'A and projected on the X-Y plane. The rot.'vtcd Y,Z coord ins^tes are 

dotarninc-d by multiplying the unrotated Y-coordinatc of the N points 

by f(A) Cos A and f(A) Sin A A respect ivsly- so dotormining another 
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function of the fanily involves detorning f(A), Sin AiA, Cos A A and 
IT+2 nultiplics-tions. If laulb iplications take tino t2 ihon lassuning 
tho-t Sin /ii A and Cos ^ A all tabs tine oqua,! to il(x) i.o., t1 j 

the detennination of anothex function of the fanilj' xequixos t2(N+2)+5t'1 
tino. The sane tine is xequixed "by all the- succossivo functions, 

Henco the total tine in Step 2 iss 

(t2(N-f2) + 5tl) * 560 /AA, let M = 560/^ A 
so the tirac in Step 2 is M(t2(N+2) + 5fcl ) 

Stop 3 loquixGS (2b 2 +3t1 ) * BMIJi/ hB ,93, , lot IIP = MM/ S © 

thon it is equal to l£P(2t2 + 3 tl). 

The fouxth step tates tine = (t1 •f-2-fc2*'l£P)*(Xt.iiX-XI;HK)/AX 
lot IIP = (Xt'IAX - X[iIir)/-^iX, then the tine equals to 
nP(t1 + 2 t 2 * IIP) 

Step 1 and 2 detexnino the fixst fanily of functions and the total 
nunhex of points on thon is and the tine xequixed is Iltl -1- M(t2(N+2)-+-3f1 
Stop 3 4 detexmine the 2nd fanily of functions vfitb tota-1 nunhex of 

points on then equal to MP^ilP and the tirae xequixod is ICP(2t2 + Jtl ) + 
lIP(t1 + 2t2 * lip). Masking loutine v;ill take appxoxinrtely 2©(M*M.!ip*IlP) 
tine (See Chaptsx 2) so the total tine xoquiied is; 

Ht1 -f M(t2(ll-f2) +3tl) + Ji!P(2t2 -t- 5t1 ) 4- iTP(t1 -fktP * Mp) + 29(n*lI-HvlP*KP) 
ox 

t2 * M ET + 2t2 MP * BP ■+• 29(M*B-<-I'ilP*‘BP) + fixst oxdox texms in 

M,B, MP and BP 

as M,B, MP, BP axe quite laxge so xetaining only the second oxdex toxns 
in then, v/e see that it is pxopoxtional to M*l'T-hI€P * BP, i.s., the total 
ninnhex of points on the suxface. 
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Inteif'cronce of Rotation SmciSaces 

Figtires 5 and 6 aie inteiference of a oylindex wibb a sphsie and 
hype 2 ;boloid of two sheet lespoctively. ‘fhe pxognani foi: singlja notation 
sircfacs is used with only a slight nodif'ication foi: t'wo 02 nore intea:- 
feiing notation sunfaces with a coanon line of rotation. In multiple 
rotation surfaces program, two masking arrays each are maintained. Wo 
take one surface at a tioo, both the families of this surf ace aie dravTO 
uding different masking arrays. This is because no member of one family 
of functions masks any part of any raember of the another family belong- 
ing to the same surface. hence the same masking arrays caiinot be 

used for both the families, ifter plotting one surface the corresponding 
tvro masking arrays (two mask low and two mask high) are combined to 
give the net masking effect of the surface plotted. Two copies of this 
resultant masking arrays are created. This would nov;’ be the base for 
plotting 2nd surface. By the very property of mamking arrays that 
portion of the second surface would be eliminatod which lies behind 
the first. So what we get essentially is projection of surfacos on 
view plane, with elimination of that portion of thera •uhich would not 
be visible to an observer at a particular angle of view. 

3-2. HAT.P TnWR RBPHBSBIlTiJTIOIi OF STJEFIOB OF THE. BOBU Z = f (x,y).. 

Half tone pictures are shedded graphic displays. The general 
technique used in half tone display is to divide the rectangle in which 
the surface is defined » into nxn squares. Depending on the local 
property of the surface in ecpb square a proper density of points 
are chosen and projected on the view plane. Size of the roctangle in 
these algorithms do not really depend whether the surface is constanffc 
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in ceitain region or fast varying, it is saim in "both the oases* 

The algorithia we have developed tries to see how the function f is 
varying, Virith X and T and depending on that a proper rectangle siae is 
chosen* The idea is to miniaise the nunher of rectangles in the region 
where the surface is defined, which would enable to reduce the time 
required to draw the surface. Once we have divided the region in 
rectangles, w.e can have various ways of shedding it. ITe have used 
the easiest way, to put a point in the centre of the reotangde. In 
the region where the function f(x,y) is slow varying we get a low 
density of points and v/here it is fast varying the points are dense, at 
the edges where the derivative, in the direction of observation is in- 
finite we get a dense cluster of points naofcing the boundary of the 
surface . 

Choosing the size of rectangles 

Suppose that the surface Z = f(x,y) is defined in the rectangle 
X-<C = X1, ^ fl* Derivatives of Z with respect 

to X and y are calculated at (x^,yj^), and the side of the square at 
th^-s point is chosen to be g/(ABS( ’J z/'S x) + XBS( "3 z/^y)) where c is 
a constant. Vfe have taken it to be 0.015 our prograa. Next the 
above procedure is repeated at (x^ + & x, y^), where 

5x = c/(iiBS( 2z / 3 x) „ + z/'3 y) ) 

o’'^o ■^o’^o 

so on until we raaobs^l • This way ire have got one row of non- 
intarese oting squares, irhose sizes are invarasedy proportional to the 
sum of absolute value of derivatives of Z with respect to x and y. The 
same procedure is repeated starting fuon 3^3. y equal bo the height 
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alxeajy aohioved atthis X-value j until the v^hols region is covered "by 
lectangles • This scheme as such leads to the piohlem cjf oveilappingj 
one rectangle inteiiferrnig with its neighbours (See Figure 7“-*^)? and also 
some space being daft uncovered between rectangles. This is avoided by 
reducing the length of the interfering rectangle to clear the overlapping 
(Figure 7~13)* It solves the interference problem but at the cost of 
X-dependencG of the length of the rectangles. Only the bottom row 
of rectangles have length depending on the derivations of the function 
f at that point, but the rectangles of the higher roYfS have length equal 
to that of the rectangle immediately below, unless the length of this 
rectangle is less than the rectangle beloYv. So virhile we traverse from 
bottom to top the number of the rectangles in a rov; will either remain 
the same or increase but in no case will it decrease, even if the deri~ 
vatives of the function decreases (see the Figure the rectangles 

ara) put one above the other, their heights vary but length remains 
the sane). To allow a certain degree of freedom in choosing the length 
of the rectangles vre have to allovY for certain amount of overlapping. 

The above scheme for interference has been modified to allow that the' 
common area botvi/een overlapping rectangles divided by the area of the 
overlapped rectangle should not exceed 0.15* If it does exceed 0.15? 
the length of the rectangle is reduced to bring it doTO to 0.15* Itie 
effeot of this row schema is seen in Figure 7”6> length of the rectangles 
can be seen to vary in the X-direction also vYith the derivatives. For 
complete algorithm see the program in Appendix II. 

Using this algorithm we have draivn sphere, paraboloid aaid hyper- 
boloid of tTiTO seats, see Figures 8, 9 and 10. 

The ticae complexity of our algorithm is dependent on the complexity 


of the surface . 



4* COITCLUSIOI aod SUGGBSTIOiJ POR FURTHER liiEPROlMvIEFrS 


Tht= peifornance of cross hatched surface algorithr. has heGn more 
or loss what had hesn intended. The a-verage GPU ticie lor about 2500 
points on the surfs.ee has been approxinately 8 SGconds-. independent 
of the surface v;hereas Watkin's algorithn takes about 19 seconds with 
2000 points for siniple surface like Z = stnx* Sin Y / XY* For conpli- 
cated surfaces like our Figures 5 4 it will take much more time, 

whereas using our algorithm they are drawn in only 6.57 and 6.5 seconds 
respecti-vely . 

The algorithm can be extended to accept multivalued H functions. 

¥g had restricted H to be single valued (see Section 5*1 )• This would 
require to modify the masking algorithm, because the present masking 
program keeps at the most two points for one pen position (they oorres- 
•pond to mask high and mask low), more than that for a particular pen 
position Yfould bo eliminated. So far multivalued functions only two of 
its single valued ranches would be kept rest of the portion would be 
eliminated even though it may be visible. 

Secondly the family of surfaces Yz-hioh can be put in the form 

= H^(x)g^(ilRCT.'ll( 2 /y)) may be extended to include those functions 
which do not belong to this family but their iavorsoe belong to it. 

For example the surface Z = 1 /x— 1 /Y cannot be put in the above form 
but i/z = be longs to the family of functions acceptable to our 

algorithm. To accept this extended family the algorithm r 6 qu,ir 0 S only a 
minor modification. All the steps are carried out exaculy as before, 
only at the time of actually drawing the surface Z-valucs are inverted* 
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j-’he half tonod surface display algoiithn h.as noi: xoally been 
a success, as the shaddsd surface does not give a good 5“!’ view of 
the surface* This algorithn is not supported hy naslciug and rotation* 
Consequently m can have only one -vievT of the surface, that is observer 
on Z-axia at infinite (even the cross hatched surface display for this 
cngln of viev^- for paraboloid and hypsrbolid docs not very much reveal 
the properties of them). Secondly for shadding v/e did not have any 
hardware support, so the dinension and intensity of the points could 
not be controlled only the density could be varied, Virithout Y/hich it is 
difficult to get a good 3"^ view of the surface. 

Possibly a better algorithm for getting shadded surfaces could be 
to choose points on the surface depending on the curvature (and not 
derivative) of the surface in its neighbourhood and then properly 
projecting it on the view plane* 
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